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I. 


INTRODUCTION 


This  program  is  aimed  at  the  development  of  radiometric 
techniques  for  the  remote  measurement  of  atmospheric  temperature 
and  humidity  profiles.  The  generated  profiles  are  expected  to 
be  instrumental  in  the  assessment  of  apparent  errors  in  target 
positions  resulting  from  refraction  index  gradients  along  an 
atmospheric  propagation  path. 

Remote  altitude  profile  measurements  of  temperature  and 
concentrations  of  gas  constituents  in  the  atmosphere  have  been 
previously  implemented  by  radiometric  measurement  of  the  up- 
welling  radiance  using  aircraft,  satellite,  and  balloon  platforms. 
(References  1,  2 and  3)-  Unfortunately  the  altitude  resolution 
obtainable  with  up-looking  radiometers  using  similar  techniques 
is  inherently  poor  due  to  the  erratic  response  of  kernel  functions 
associated  with  down-welling  radiation.  Therefore,  the  measured 
information  for  an  up-looking  radiometer  lacks  sufficient  detail 
to  retrieve  properties  of  the  lower  atmosphere  where  sharper 
gradients  of  temperature  and  humidity  often  occur  as  a function 
of  diurnal  and  other  atmospheric  conditions. 

An  alternate  technique  for  obtaining  more  accurate  tempera- 
ture and  humidity  profiles  with  a grazing-angle  up-looking  radiom- 
eter has  been  investigated.  The  profiling  technique  is  based  on 
information  retrievals  which  are  obtained  from  difference-kernel 
functions  generated  from  pairs  of  high  spectral  resolution  radio- 
metric  channels.  Spatially  selective  information  can  be  generated 
from  such  functions  with  the  required  accuracy  for  retrieving 
temperature  gradients  at  the  lower  altitude  atmospheric  levels. 
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A retrieval  of  an  assumed  double  inverted  temperature  profile 
in  the  lower  atmosphere,  using  three  differential  kernel  information 
channels  of  COg  signature  emission  in  the  13  pm  region,  has  been 
simulated.  The  retrieval  is  based  on  iterative  solutions  of  the 
radiative  transfer  equation  and  requires  no  apriori  statistical 
knowledge  of  the  atmospheric  temperature  profile.  Temperature 
inversions  at  375  and  750  meters,  assumed  for  the  selected  atmos- 
pheric model,  were  reproduced  in  the  simulated  retrieval.  However, 
the  retrieved  profile  lacked  detail  at  the  immediate  altitude  above 
ground  level  due  to  the  excessive  spatial  widths  of  the  difference- 
kernel  functions. 

A four  information  channel  retrieval,  based  on  higher  selec- 
tive functions , will  be  investigated  during  the  next  phase  of  this 
program. 

II.  GENERAL  SCOPE 

In  general,  variations  in  the  atmospheric  refraction  index 
may  be  classified  as: 

• Small-scale 

* Large-scale 

Small  scale  variations  are  associated  with  clear  air  turbulence 
phenomena  in  the  lower  atmosphere  and  are  manifested  in  random  fluc- 
tuations of  amplitude  and  phase  of  energy  transferred  through  the 
medium,  with  a time  constant  in  the  millisecond  range.  The  large 
scale  refraction  index  profile;  of  interest  in  this  study,  is  a 
relatively  long  time  property  of  an  atmospheric  path  which  is  a 
function  of  average  values  of  local  pressure,  temperature  and 
humidity  along  the  path. 
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The  refractive  index  n of  the  atmosphere  is  conveniently 
described  in  terms  of  the  ref ractivity  N,  in  parts  per  millions, 
by  the  expression, 


n = 1 + 10  n N 


or  N = (n-1)  10c 


The  refractivity  at  optical  wavelengths  is  related  to  local 
temperature  and  pressure  (Reference  4)  by: 

N ^79  | (2) 

where:  P = pressure  in  millibars 

T = temperature  in  degrees  kelvin 
The  downward  bending  curvatuve  of  an  optical  ray  at  a given 
altitude  level  is  given  by  the  gradient  of  the  refractivity: 

„ dN  79  dP  , 79P  dT 


For  radio  and  microwave  frequency  propagation,  an  additional  term 
due  to  humidity  gradient  enters  into  the  above  curvature  equation. 

Absolute  values  of  local  temperature  and  pressure  as  a func- 
tion of  altitude  entering  in  equation  (7),  may  be  estimated  within 
a small  error  from  reported  statistical  data  for  given  geographic 
area  and  measured  ground  (or  ocean)  temperature.  Similarly,  close 
estimates  on  pressure  vs  altitude  gradients  are  obtained  from  univer- 
sal curves.  The  situation  however  is  quite  different  for  the  case 
of  temperature  gradients.  Available  statistical  data  on  temperature 
profiles  provide  good  estimates  on  gradients  at  relatively  high 
altitudes.  At  the  lower  altitudes  however,  changes  in  the  gradient 
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directions  (i.e.,  inversions)  occur-  frequently  as  a result  of  many 
atmospheric  variables,  (e.g.,  position  of  the  sun,  overcase  condi- 
tion, surface  and  ambient  temperature,  wind  velocity  and  direction), 
which  result  in  substantial  error  in  assessing  the  true  target  position. 

Measured  temperature  gradients  at  lower  altitude  using  radio- 
sonde techniques  are  given  in  Figure  1.  The  double  inversion  in 
the  gradient  (superadiabatic ) is  typical  for  mid-day  condition 
whereby  the  surface  temperature  is  higher  than  the  surrounding 
ambient  temperature.  Positive  inversion  at  ground  level  (not  shown 
in  Figure  1),  on  the  other  hand  typically  occurs  at  night.  The  neutral 
gradient  (Figure  1)  is  normally  typical  for  early  morning  and  later 
afternoon  conditions.  It  is  clear  that  the  uncertainty  in  the 
changing  gradient  and  associated  inversions  in  vicinity  of  ground 
or  ocean  surface  is  the  largest  source  of  potential  error  in  esti- 
mating the  local  curvatures  in  ray  bending  for  grazing-angle  up- 
looking  pointing  or  direction  finding  detection  systems. 

'This  study  is  aimed  at  establishing  the  feasibility  of  deter- 
mining the  temperature  and  humidity  gradient  profiles  at  the  lower 
altitude  levels  by  means  of  remote  radiometric  measurements.  Using 
this  measurement  information  in  conjunction  with  statistical  'measured 
data  on  local  atmospheric  conditions  (temperature,  pressure  and 
humidity),  assessing  target  error  positions  in  real  time  appears 
to  be  feasible. 

III.  BASIC  PHENOMENOLOGY  IN  RADIOMETRIC  PROFILING 

The  average  apparent  radiance  sensed  by  a radiometer  viewing 
the  atmosphere  is  governed  by  the  steady-state  atmospheric  radiation 
transfer  equation  as  illustrated  schematically  in  Figure  2. 


For  a down- looking  radiometer  (Figure  2A),  the  apparent 


radiance  is  given  by: 


N 


A 


3i(T) 


"JT. 
O 1 

V 

o J 


dy 


where : 


(M 


BQ,  B^.  = Planck's  radiances  at  ground  and  ith  altitude-- 
level,  respectively 

t , t.  = atmospheric  transmission  coefficients  in  the 
total  propagation  path  from  ground  level  to 
receiver,  and  ith  altitude  level  to  receiver, 
respectively 

y = normalized  spatial  coordinate 

The  first  term  in  the  right-hand  side  of  equation  (4)  is  due  to  the 
earth  (or  ocean)  radiance  as  attenuated  in  its  total  propagation 
path.  The  second  term  represents  re-emission  of  radiance  in  the 
various  atmospheric  layers  within  the  total  propagation  path. 

A similar  transfer  equation  may  be  obtained  for  an  up- looking 
radiometric  system  (Figure  2T. 


NA  - J Bt<T>  — (5) 

yt 

where  the  sensed  radiance  is  due  to  atmospheric  emission  only. 

One  method  of  radiometric  profiling  consists  of  implementing 
several  spectral  sensor-channels , each  selectively  accepting  emitted 
radiation  from  limited  spatial  regions.  Such  an  ideal  up-looking 
profile  radiometric  system  is  depicted  in  Figure  3.  As  shown,  it 
consists  of  a multitude  of  sensor  channels,  each  having  rectangular 


Figure  3.  Ideal  Profiling  - Radiometer  Concept 


• •• 


narrow  response  characteristics  sensitive  to  radiation  from  dis- 
crete altitude  levels.  In  practice,  the  responses  of  the  channels 
are  broadened  by  natural  characteristics  in  the  signature  profiles 
of  the  emitting  gases  as  well  as  by  their  relative  distribution  with 
altitude  levels.  Since  the  response  widths  of  the  actual  channel 
are  finite,  they  must  be  contiguous  so  as  to  encompass  the  total 
profiled  range  of  interest. 

The  basic  parameters  associated  with  radiometric  profiling 
are  defined  in  Figure  4,  where  a down-looking  radiometer  is  view- 
ing an  "atmospheric  layer"  within  its  field-of- view.  The  pertin- 
ent radiometric  parameters  are: 

B^(T)  = Planck's  radiance 
Wj_(p)  = concentration  of  active  gas 
£(w,P)  = gas  emissivity 
T = temperature 
p = pressure 

It  should  be  noted  that  though  gas  concentration  and  emissivity  are 
also  functions  of  temperature,  this  functional  dependency  is  small 
as  compared  to  that  due  to  pressure,  and  hence,  it  is  neglected  here. 
The  apparent  radiance  emanating  from  the  layer  as  sensed  by  the 
radiometer  is  shown  in  Figure  4 and  given  by: 

% = Bi(T)  €l-zi  (6) 

where  x^  is  the  transmission  coefficient  of  the  total  path  from  ith 
level  to  receiver. 
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The  apparent  emission  relationship  in  equation  (6)  may  be 
expressed  in  a somewhat  alternate  form.  Referring  to  Figure  4, 
assume  a radiance  source  B (T)  alternately  located  at  points  A 
and  B,  below  and  above  the  atmospheric  layer  considered,  respec- 
tively. The  apparent  absorption  within  the  layer  is  then  given  by: 


From  Kirchof's  law  of  absorption-emission  equivalency,  the  coeffic- 
ient in  equation  (6)  is  equivalent  to  the  transmissivity  difference 
in  equation  (7)  which  may  also  be  expressed  in  a differential  form, 
namely : 

aT. 

iTi  = T1  " t2  — JJ 

STi 

The  spatial  derivatives  — -r — -,  referred  to  as  kernels  or  weighting 

<j  y 

functions,  are  the  linear  coefficient  of  the  temperature  dependent 
Planck's  radiance  in  equations  (4)  and  (5)»  Hence,  the  spatial 
selectivity  of  a radiometer  is  determined  from  response  character- 
istics of  the  kernels  as  a function  of  altitude. 

The  difference  in  the  nature  of  response  of  the  kernels  for 
up-looking  and  down-looking  systems  may  be  visualized  by  referring 
to  Figure  5-  For  a down- looking  radiometer  (Figure  rjA ) , the  con- 
centration in  the  active  gas  and  hence  the  emissivity  is  decreasing 
along  the  direction  of  radiation.  The  transmission  coefficient 
increaases  along  the  direction  of  the  propagation  pat,}:.  Hence, 


dy 


(8) 


peaking  of  the  weighting  function  along  the  propagation  path  will 
occur.  For  an  up-looking  radiometer,  no  peaking  will  occur  since 


both  quantities  increase  along  the  direction  of  radiation  as  shown 
in  Figure  5B.  Typical  monochromatic  weighting  functions  as  a func- 
tion of  altitude  are  shown  in  Figure  6 for  a down-looking  system, 
and  Figure  7 for  an  up-looking  system.  As  may  be  seen  from  Figure  7, 
the  spatial  resolution  for  an  up-looking  radiometer  is  poor  due  to 
the  exponential  nature  of  the  weighting  functions. 

IV.  DIFFERENTIAL  KERNEL  APPROACH 

In  view  of  the  above  discussion  it  is  clear  that  weighting 
functions  generated  in  channels  of  an  up-looking  radiometer,  which 
are  solely  based  on  spectral  variations  in  the  signature  strengths 
of  the  active  gases,  suffer  from  inherently  poor  spatial  selectivity. 

Hence,  additional  properties  in  the  signature  parameters  must  be 
utilized  to  generate  selective  kernel  functions.  A possible  approach 
is  the  utilization  of  the  inherently  steep  slope  characteristics 
associated  with  the  individual  rotational  emission  lines  of  atmos-  1 

pheric  gases  as  an  additional  degree  of  freedom.  This  requires 
sensing  radiation  from  sections  of  a rotational  emission  line  by 
means  of  a high  spectral  resolution  radiometer. 

Figure  8 illustrates  a rotation  line  shape  with  Lorentzian 
attenuation  profile  which  well  represents  pressure  regions  in  the 
vicinity  of  ground  levels  of  the  atmosphere.  fhe  absorption  as  a 
function  of  wave  number  is  given  by: 

r 


ALTITUDE  (KM) 


ABSORPTION  COEFFICIENT 


Figure  8.  Lorentz  Line  Shape  Profile 


where 


S = line  intensity 

T - half-width  at  half  maximum 
v>v0  = variable  and  line  center  wave  numbers,  respectively 

The  altitude  dependency  of  the  absorption  coefficient  is  due  to  the 

pressure  variation  in  line  wridth  and  is  given  by: 

r = rQp  (l0) 

where 

p = = normalized  pressure 

o 

PQ  = pressure  at  ground  level 

y = line  width  at  ground  level 
' o ° 

The  altitude  (pressure)  dependent  transmission  coefficient  for 
monochromatic  radiation  (which  is  a good  approximation  for  high 
resolution  radiometers  with  channel  spectral  widths  of  the  order 
or  less  than  1 percent  of  the  absorption  line  width)  is  given  by: 


w(p ) = wQp 

where 

w ( p ) = gas  concentration 

wq  = concentration  at  ground  level 
a(p)  = absorption  coefficient 
0 = zenith  angle 


h = altitude 


(11) 


(12) 
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Inserting  equations  (9),  (10),  and  (12)  into  equation  (11), 
and  differentiating  with  respect  to  scale-height  parameter  y,  the 
following  expression  for  the  weighting  function  is  obtained: 


T 

y 


k sec  0 
'(l+C2)*''2 


2 2 
C +p 


(15) 


where 


Sw 


k = 


7T 


r o 

V- V 


C = 


Parameters  independent 
of  pressure 


p = e 


-y 


; ^y  x 


h ( km ) 


5 


pressure  scale  height 
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As  an  example  consider  the  functional  relationship  of  the 
weighting  function  with  pressure  for  extreme  spectral  positions  of 
sensor  channels  at  two  rotation  lines  of  various  signature  intensities. 


A.  For  the  first  channel  located  at  line  center: 


C = 0 


_£_t 

py 


(14) 


B.  For  the  second  channel  located  at  the  far  wings  of  the  line: 


C >>  p 


py 


(15) 
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The  ait ference-kernel  generated  from  the  absorption  lines  with  cor- 
responding channel  positions  (equations  11  and  12)  is  given  by: 

= Apk  - Bp2  fl6) 


with  a peak  occurring  at  pressure  level  altitudes : 


'peak 


r 2 Bi 

A J 


where : 


k sec  Q 
(1+C2)k/2 


k.  sec  9 

(i+c1")ki/" 

Thus  peaked  weighting  functions  may  be  generated  from  difference 
signals  at  two  spectral  line  positions,  with  a spatial  selectivity 
being  a function  of  the  radiative  signature  intensities  of  the  two 
lines.  It  should  be  pointed  out,  however,  that  the  difference 
kernel  obtained  from  channels  at  extreme  spectral  positions  are 
rather  broad  due  to  the  shallow  gradient  at  the  far  wing  of  the 
line  and  hence  are  inadequate  for  profiling  the  boundary  layer 
with  the  required  high  spectral  resolution. 

Referring  to  Figure  8,  it  may  be  observed  that  a steep  slope 
in  the  absorption  line  occurs  in  the  vicinity  of  its  half  intensity 
position  and  hence  its  sensitivity  to  pressure  variation  will  be  large 
in  this  region.  This  suggests  that  narrower  response  difference  kernels 
may  be  generated  for  channel  positions  in  the  vicinity  of  its  steepest 
slope.  This  configuration  is  shown  in  Figure  Q. 
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Figure  9* 


Schematic  2-Channel  Arrangement  for  Generating  Difference 

Kernels 
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GENERATION  OF  DIFFERENCE  KERNEL  FUNCTIONS 


The  COg  emission  band  at  13  pm  was  considered  for  generating 
differential  kernel  functions  due  to  its  high  signature  intensity 
characteristics  and  the  feasibility  of  high  spectral  resolution 
instrumentation  in  the  band  (e.g.,  heterodyne  radiometric  techniques 
using  tunable  diode  laser  local  oscillators).  The  technique  however 
can  be  adapted  at  other  spectral  bands  in  the  infrared  or  millimeter 
waves  where  discrete  signature  spectra  due  to  uniformly  mixed  atmos- 
pheric gases  exist. 

The  difference  in  kernels  for  two  channels  associated  with  two 
discrete  lines, is  given  (see  equation  (13))  by: 


* * Hyb 


-ii) 
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sec  0 > K 
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L i+c. 
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K2 

p 

. 1+C2 

c22+<- 

(18) 


where : 


Kf,Kg  relate  to  signature  strength  of  the  two  lines  as 
defined  in  equation  (13) 

C,,Cp  are  normalized  spectral  positions  of  the  two  line 
channels 

Thus,  given  two  spectral  lines  of  signature  strengths  and  K , the 
spectral  positions  of  the  channels  0^  and  C0  which  result  in  peaking 
the  difference  kernel  at  a given  altitude,  are  found  by  minimizing 
equation  (l8)  at  a pressure  corresponding  to  the  altitude  level, 
namely : 
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Such  computations  were  performed  by 


and  solving  for  C-^  and  C0 . 
solving  the  transcendental  equation  (l8)  for  various  combinations 
of  line  pairs  and  channel  spectral  positions,  so  as  to  obtain  highly 
selective  and  contiguous  difference  kernel  channels  at  the  lower 
atmospheric  levels. 

Figures  10  and  11  show  difference-weighting  functions  ob- 
tained for  vertical,  and  30-degree  grazing-angle  up-looking  sys- 
tems, respectively.  As  shown,  narrow-width  difference  kernels 
result  at  contiguous  regions  immediately  above  ground  level.  The 
dotted  lines  shown  in  the  figures  represent  undesirable  reversal 
in  sign  (foldovers)  in  the  difference  functions  which  must  be  dealt 
with  in  the  retrieval  analysis  so  as  to  minimize  this  potential  erro 
A simulated  temperature  retrieval  based  on  the  difference  kernel  was 
carried  out  and  is  discussed  in  Section  IV. 

VI.  COMPUTER  SIMULATED  RETRIEVAL  OF  SUPERADI ABATI C TEMPERATURE 

Pit  -TIE 

Retrievals  of  a typical  superadiabatic  lapse  temperature 
profile  were  computer  simulated  using  the  three  lower  difference- 
kernel  information  channels  in  Figure  11.  The  algorithm  adapted 
is  based  on  iterative  solutions  of  the  radiative  transfer  equations 
and  require  no  apriori  statistical  knowledge  on  atmospheric  param- 
eters. The  method  is  based  on  the  up-welling  radiation  retrieval 
proposed  by  W.  L.  Smith  (Reference  5)  which  was  modified  here  for 
a down-welling  situation.  The  basic  outline  of  the  iterative 
computation  procedure  is  as  follows. 
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A temperature  altitude  profile  is  first  assumed,  and  desig- 
nated as  the  "true"  model  atmosphere.  The  radiance  corresponding 
to  the  "true"  profile  at  each  spectral  channel  is  then  computed 
from  equation  (5)  by  the  method  of  digitization  of  the  atmosphere 
in  discrete  isobaric  incremental  regions.  The  computed  radiance 
levels  at  the  channels  are  designated  as  the  "measured"  values. 

An  arbitrary  temperature  profile  designated  as  a "guessed"  model 
is  then  assumed  serving  as  a starting  point  in  the  iterative  comp- 
utational process.  Hie  radiances  corresponding  to  the  "guessed" 
model  for  all  spectral  channels  as  well  as  radiance  ratios  of 
"measured"  to  "guessed"  values  are  initially  computed. 

As  a first  step  in  the  iteration,  "guessed"  model  radiances 
corresponding  to  each  incremental  altitude  layer  are  updated  in 
direct  proportion  to  the  computed  radiance  ratios  of  "measured" 
to  "guessed"  values.  A first  iterated  model  is  then  generated  by 
averaging  the  temperatures  corresponding  to  radiances  of  the  spectral 
channels,  each  weighted  by  their  characteristic  kernel  function  at 
the  given  altitude.  A relaxation  method  of  iteration  is  adapted 
wherebv  spectral  radiances  at  incremental  atmospheric  layers  are 
continuously  updated,  proportionately  to  the  ratio  of  the  "measured" 
to  computed  values.  New  consecutive  models  are  generated,  based  on 
averaged  weighted  temperature  parameters  from  differential  kernel 
functions  at  the  given  altitudes.  The  process  continues  until  the 
difference  between  "measured"  and  computed  spectral  radiances  are 
below  a predetermined  threshold  value. 
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A simulated  temperature  retrieval  for  an  assumed  double 
inverted  temperature  profile,  based  on  three  information  channels 
(Channels  1,  2,  and  3 in  Figure  11)  is  given  in  Figure  12.  Curve 
A in  the  figure  is  the  assumed  "true"  atmospheric  profile,  the 
dashed  curve  is  the  arbitrary  assumed  initial  model  in  the  itera- 
tion, and  the  retrieved  model  profile  is  as  shown  in  curve  C. 

The  two  inversions  in  the  assumed  model  at  37^  m and  750  m alt 
tude  are  reproduced  in  the  retrieval.  However,  the  fidelity  of  the 
retrieved  profile  in  the  lower  375  meters  above  ground  level  was  in- 
sufficiently inaccurate  due  to  excessive  spatial  widths  in  the 
kernel  functions  of  the  three  channels.  A more  accurate  tempera- 
ture retrieval  implementing  four  information  channels  with  narrower 
spatial  functions  will  be  investigated  during  the  next  phase  of  the 
program. 

VII.  INFRARED  HIGH  SPECTRAL  RESOLUTION  RADIOMETRY 

The  required  high  spectral  resolution  of  10~  to  10-J  cm"^ 
for  sensing  radiation  from  portions  of  a rotational  emission  line 
necessitates  the  utilization  of  state-of-the-art  spectrometer  tech- 
niques. Present  day  conventional  spectrometers  do  not  provide  the 
required  resolution.  However,  a heterodyne  approach  may  be  employed 
wherein  a tunable  laser  (diode  local  oscillator)  illuminates  a photo- 
diode simultaneously  with  the  incoming  source  irradiance.  The  result 
ing  IF  signal  is  a high  spectral  resolution  slice  of  the  infrared 
radiation.  This  technique  has  been  successfully  employed  in  a number 


of  equipment  and  instrumentation  developments  (References  2 and  3)- 


A dual-channel  Infrared  heterodyne  radiometer  (IHR),  as 
shown  in  Figure  13  (Reference  6),  is  currently  being  used  to 
assess  the  concentration  and  vertical  distribution  of  atmospheric 
ozone  and  ammonia  by  observing  the  up-welling  radiance  from  differ- 
ent parts  of  the  selected  gas  constituent  line.  The  heterodyne 
radiometer  techniques  employed  to  measure  gas  concentrations  can 
also  be  used  to  retrieve  the  thermal  profile. 

VIII.  SUMMARY 

A new  method  of  remote  profiling  temperature  and  humidity 
gradients  in  the  lower  atmosphere  by  means  of  up-looking  grazing- 
angle  radiometer  has  been  investigated.  The  technique  consists  of 
generating  difference-kernel  functions  in  several  high  spectral 
resolution  channels.  These  functions  are  shown  to  have  inherent 
high  spatial  resolution  in  the  lower  regions  of  the  atmosphere. 

The  results  of  a simulated  retrieval  of  temperature  profile  with 
an  assumed  double  inversion  below  an  altitude  of  1 km,  based  on 
three  information  channels,  is  presented. 
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